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\Q ', Abstract 

I analyze the one-dimensional, cubic Schrodinger equation, with nonlinearity con- 
structed from the current density, rather than, as is usual, from the charge density. A 
soliton solution is found, where the soliton moves only in one direction. Relation to 
higher-dimensional Chern-Simons theory is indicated. The theory is quantized and 
results for the two-body quantum problem agree at weak coupling with those coming 
from a semiclassical quantization of the soliton. 



1 Introduction 



The Schrodinger equation in one spatial dimension with cubic nonlinearity 

d -h 2 d 2 

ih—ip(t,x) = - — -^ip(t,x) - gp(t,x)ip(t,x) (1) 

g > 0, p = ip*tp . 

plays a cycle of interrelated roles in mathematical physics. (Here * signifies complex 
conjugation.) 

Viewed as a nonlinear partial differential equation for the function i/j, it possesses the 
famous soliton solution 

Mt,x) = ±^(M » « (2) 

y/gm cosh a [x — vt) 

a 2 = (m 2 v 2 /h 2 )(l - 2u/v) . 
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Moreover, Eq. (|l|) is completely integrable and multisoliton solutions can also be explicitly 
constructed. Note that the soliton described by Eq. (||) moves with (group) velocity v, 
which is unrestricted as to sign and magnitude. (The phase velocity u must be less than 
v/2.) In particular the soliton can be brought to rest: the theory is Galileo-invariant, so 
that any solution to (|l|) can be mapped into another solution by the Galileo boost 

x ► x — Vt 

i>(t, x) ► e $ mV <- x -* v H(t, x-Vt) . (3) 

[Galileo transformations for the dynamical system (|l]) are accompanied by a phase change 
of ijj - a so-called "1-cocycle" .] Applying the transformation (@) to (g) with V = —v 
results in 

i> a (t, X) ► ±e l 2m —— 4) 

yjgm cosh ax 

which also solves ([!]) and describes a soliton at rest. 

Alternatively one may view (|l]) as a quantal Heisenberg equation of motion for the 
quantum field operator ip(t, x), which is taken to satisfy the commutation relation 

[ip(t,x 1 ),ip*(t,x 2 )] = 8(xi - x 2 ) ■ (5) 

(Now * denotes Hermitian conjugation.) Since the dynamics (jl]) conserves the number 
operator J dx p, the quantum Hilbert space can be decomposed according to the (integer) 
eigenvalues iV of J dx p, and one finds that the A^-body wave function 

ip N (t;xi,. . .,x N ) = -^==(0\i/;(t,xi) ■ ■ ■ tf>(t, xn)\N) (6) 

satisfies an A-body Schrodinger equation with two-body, pairwise attractive (^-function 
interactions. The bound-state spectrum can be found explicitly and the energy eigenvalue 
is determined as 

E N = ^(N 3 -N). (7) 

The cycle of formulas closes: the static solution (|j) can be quantized semiclassically 
and the bound energy spectrum obtained in this way coincides with (]7j).[]l[] Specifically, 
the classical energy of (||) is A 3 , where A is the value of J dx p evaluated on (^). 
The semiclassical (first quantum) correction diminishes A 3 by A (so that the bound state 
energy vanishes at A = 1). (Also, the semiclassical phase shift of the scattering problem 
reproduces the exact quantal formula for weak coupling. |J) 

In this article, I shall review some properties of another nonlinear Schrodinger equation, 
which (in a transformed formulation) has recently arisen in discussions of one-dimensional 
condensed matter systems (quantum wires, Hall edge states) || and which possesses a 
remarkable soliton solution, in that the soliton is chiral, i.e., it can move only in one 
direction.^] 

2 The Equation and Its Soliton Solution 

The equation that we consider differs from ([l]) in that the nonlinearity involves the current 
density j 

j = -Im(^) (8) 

m \ ox ' 
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rather than charge density p: 

d h 2 d 2 

ih—ip(t, x) = --^—-^-o^it, x ) - K*>j(t, x)ip(t, x) . (9) 
ot 2m ox 

Note that j can be substituted for p only in one spatial dimension, where the "vector" 
nature of j is absent. The current and charge densities are linked by the continuity 
equation, which is valid both for (|l|) and @. 

d d 

ai» + TJ = a - (10) 

When if) is decomposed into modulus and phase 

$ = y[p&* (11) 

one sees that 

'■ = £4' < 12 > 

and ([|) takes the form of the usual nonlinear Schrodinger equation (||), 

d -h 2 d 2 

ih— ip(t, x) = 2m dx 2 ^' ^ ~ 9<(t ' X ^ P ^' ^ ( 13 "> 

but with coupling strength 

g(t,x)=h 2 ~9(t,x) (14) 
m ox 

modulated by the spatial variation of the phase. There is no a priori restriction on the 
sign of A, but for definiteness we shall take it to be positive, A > 0. 

To find the one-soliton solution tp s , we take the phase to be as in a propagating wave 

^ = e _ iM _ M ^_ (15) 

so that 

j = vp, v = — (16) 
m 

and our equation (g) becomes identical to ([j]), with 

g = h\v . (17) 

The soliton solution exists provided g > 0; this requires 

v > . (18) 

The soliton can only move to the right; it is chiral. (Had we taken A < 0, then the soliton 
could move only to the left.) The profile that solves @ is found from @ and (|l7| ) to be 



ijj s (t,x) = ± e i ^ x - ut \ — - (19) 

V Xmv cosh a(x — vt) 

where u = uj/k. The presence of 1/y/v in the amplitude reinforces the statement that the 
velocity must be positive; in particular, the soliton cannot be brought to rest. Evidently, 
the dynamics @ is not Galileo- invariant; the velocity of the soliton cannot be arbitrarily 
reduced. This fact will be seen clearly in the next section, where the action, Lagrangian, 
and other dynamical quantities relevant to @ are discussed. 
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3 Action Principle, Symmetries, and Constants of Motion 

Eq. (|9|) may be obtained from an action principle, which also is useful for identifying the 
Hamiltonian (energy), momentum, and other constants of motion. 
Consider the Lagrange density 



C 



a 



iW— * 

ot 2m 



The Euler-Lagrange equation reads 

Ot 



2m ^dx 



it I a A \2 h\ ., 
1 ~P) ^-—3® 



where 



(20) 

(21) 
(22) 

V va, -2-J-j (23) 

and the two are linked by the continuity equation ([n]). Next, we redefine the \P field by 

y(t, x) = e** /* dy ^V(i, s) (24) 
(the lower limit on the integral is immaterial - it affects only the phase of tp) so that 



P 

3 



— Im (&* (S- - i~p) * 

m 



X 



d 



TiX 



dt 







at' 



~^il,(t, x) - ^j(t, x)1>(t, x) . (25) 
2m ox 2 



But the integral may be evaluated with the help of ([!(]) and transferred to the right side. 
The resulting equation is just (||). 

Energy is conserved as a consequence of time-translation invariance, and its form can 
be deduced with the help of Noether's theorem, or by inspection from (|20|). Evidently the 
Hamiltonian (energy) density is 



£ 



2m 



(- 

V<9x 



.A \ 2 



2m 



<9x 



(26) 



Similarly, space translation invariance ensures momentum conservation, and the momen- 
tum density reads 



7> = ftIm(V— V) 



?iA 2 
mj + — p 



(27) 



These quantities obey continuity equations. With the energy density £ there is associated 
an energy flux T 0x 



-.Ox 



d _r d 



m 



m 



iX 



dx 2 p 



•) 



t?.a 



(28) 
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and together they satisfy 



— £ + —T 0x 
dt dx 



. (29) 
Similarly, with the momentum density V, supplemented by the momentum flux 



K 2 

rpXX 

m 
m 



, d «A 



2 n 2 d 2 



Am dx 2 



P 



— / 

dx 



tf_cP_ 
Am dx 2 



P 



their continuity equation is verified 



d m __ 
-p _| t 

dt dx 







(30) 



(31) 



Therefore, when fields decrease rapidly at spatial infinity, the energy and momentum are 
time-independent 



E 
P 



dx£ ^-E 
dt 







dxV 



dt 



-P = . 



(32) 
(33) 



Note that the momentum density V is not equal to the energy flux T 0x , because the 
theory is not Lorentz-invariant. Also it is not Galileo-invariant. This is seen in the present 
context by constructing the usual Galileo boost generator 



Prom fll0|) and (|33]) we find 

dG 
~dt 



G = tP-mJdxx P . 
P + m J dxx —j = P — m J dxj. 



(34) 



(35) 



But Eq. ( g7j ) shows that P possesses the dynamical contribution — J dx p in addition to 
the usual kinematical term m J dxj; therefore it follows that G is not conserved, 

d 



dt 



G 



~2~ 



dx p 



(36) 



but always increases in time. 

There does exist, however, another symmetry, beyond time and space translation in- 
variance and another conserved quantity beyond H and P. The additional symmetry is 
dilation invariance, corresponding to the possibility of rescaling space and time: 



ip(t,x) 



-> at, x — 
— > a 1/ V(at 



* yax, 
ax) . 



(37) 



Applied to the soliton solution (|T^) , the transformation ( |37| ) yields another soliton solution 
with group and phase velocities v, u rescaled by y/a. In order to construct the constant of 
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motion associated with dilation invariance, it is useful to "improve" the energy density and 
energy flux by the addition of terms ( "superpotentials" ) that do not affect the continuity 



equation nor the spatial integral defining the energy. Instead of fl2q) and (25), we use 



h 2 c?_ 

dx 

n 2 d 2 



improved -£ Smdx 2P ( 38 ) 



improved 8m dx^' 



Owing to (0), we still have 

d_ | d_ 

^■"-improved + ~~Q^. ± improved 

and for fields that decrease at large distances 



^improved + n ^Improved ~~ ^ (^0) 



E = J dx£ = J dxS impiaved . (41) 
The improved energy density satisfies 

^improved — 2'P (^^) 

which is the nonrelativistic criterion of scale invariance in a field theory. M It follows from 
(|42| ) that the dilation change 

D = tE-±JdxxV (43) 

is time-independent: 

^ = E-±fdxx£ i T> = E+±fdxxJ^T xx 

= J dx(£ impmved - \T**) = . (44) 

[Condition (fS) also ensures conformal invariance in a Galileo- invariant theory. || Here 
the absence of Galileo invariance requires absence of conformal invariance, as is seen from 
the Lie algebra of the corresponding generators: the bracket of P with the conformal 
generator closes on the Galileo boost generator.] 

The momentum and energy of the soliton ( |T^ ) are obtained by evaluating P and E on 
the profile 

Ronton = Mv (45) 

Ronton = \Mv 2 = ^P 2 ohton (46) 
/ A 2 \ 

M = mN(l + — iV 2 J (47) 

with 

N = J dxp=|^. (48) 

The soliton's dynamical characteristics are those of a nonrelativistic particle of mass M 
moving with the velocity v. The consequent connection between energy and momentum 
exhibited by @ and (g§), 

-^"soliton 2 ^^"soliton (^^) 
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in fact follows dilation invar iance. For the soliton fly]), V is a function x — vt, hence the 



time-independent dilation change (43) reads 

D = tE-\l dxxV(x - vt) 

= t[E- ~p) - ±fdxxV(x) (50) 

where the second equality follows from the first by a shift of integration variable. Since 



both D and the last term in (4£) are time-independent, the coefficient of t must vanish, 
and this implies (|49|). Since for our solution V(x) is an even function of x, the last integral 
also vanishes and 

Aoliton = . (51) 

4 Other Solutions 

Thus far we have considered only the attractive interaction, which binds a soliton: g > 
in (jl|) and Xv > in (17, |T1| ), With a repulsive interaction, there also exist classical 



solutions, variously known as "kinks" or "dark solitons" . U They carry infinite energy, 
and we shall not discuss them here. 

Unlike the conventional nonlinear Schrodinger equation ((!]), our chiral equation (g) 
does not appear to be completely integrable Q, and analytic expressions for multisoliton 
solutions to (||) are not available. It has been remarked, however, that if the theory is 
modified by adding to © the potential V(p) = -^-p 3 , then Eq. (|J) acquires the extra 
term ~ 3 g^ p 2 ^ and becomes an integrable nonlinear, derivative Schrodinger equation 
with nonlinearity i\^P~§z^- However, the solitons are no longer chiral.^] 

5 Relation to (2+l)-dimensional Chern— Simons Theory 



Our dynamical, chiral model with Lagrange density (|20| ) is partially related, through 
dimensional reduction, to a (2 + l)-dimensional model of nonrelativistic fields interacting 
with a U(l) gauge potential, whose kinetic term is the Chern-Simons expression. 
Consider the (2 + l)-dimensional Lagrange density 

2 2 



i=i 

When dependence on the second spatial coordinate is suppressed, and A2 is renamed ^tB, 
£(2+1) leads to a B-F gauge theory: 

Ai+D = y*#" F i»> + im * (i + * - £ 1 [tx + u -)*r - w B2p • (53) 



2k ^ \dt 'J 2m 

9 A 9 a rvM _ {+ „\ A — ( A~ A \ v — & 



Here = -^pA v - ■$£?A lt , x^ = {t,x}, A^ = (A ,A X ), and k = — k, where c is the 
velocity of light, which plays no role in the following. This is not yet our theory (^0[), be- 
cause when the nonpropagating B and A^ fields are eliminated, they decouple completely 
in the sense that the phase of \£ may be adjusted so that all interactions disappear: 



2 



dx 



2 



(54) 
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In order to make the vector potential and the B field dynamically active, thereby 
allowing the \& particles to interact, we include a kinetic term for B, which could be taken 
in the Klein-Gordon form. However, we prefer a simpler expression that describes "crural" 
Bose fields, propagating only in one direction, whose Lagrangian density is proportional 
to ±^§§ + ^^§^§-[|9[ Here v is a velocity and the consequent equation of motion arising 
from this kinetic term (without further interaction) is solved by B = B(x ^ vt) (with 
suitable boundary conditions at spatial infinity), describing propagation in one direction, 
with velocity ±v. Note that is not invariant against a Galileo transformation, which 

is a symmetry of £(1+1) and of performing a Galileo boost on with velocity 

v gives rise to vQj^-^, effectively boosting the v parameter by v. Consequently one can 
drop the v f|rf|r contribution to the kinetic B Lagrangian, thereby selecting to work in a 
global "rest frame" . Boosting a solution in this rest frame then produces a solution to the 
theory with a §f §f term. 

In view of the above, we choose the .B-kinetic Lagrange density to be 

ldBdB 

and the total Lagrange density is Cb + It is possible to remove the and B 

fields by a Hamiltonian reduction, as described in Ref. [fief ], and by a phase redefinition of 
ty. Once this is done, one is left with the Lagrange density ( p0|) 

£b + £(i + i) >£ (56) 

with ±k 2 entering as — j. The "±" sign reflects the sign arbitrariness of the kinetic 
term, and without loss of generality can be chosen so that A > 0. Also, the origin of 
Galileo noninvariance is now recognized as arising from £ B .0 



6 Quantum Theory 

The theory described by the Lagrange density (^) may be quantized. For the quantum 
Hamiltonian we take the normal ordered expression 



h 2 

H = — I dx 

2m 



fd .A x 2 
(to- 2')* 



(57) 



and posit the canonical commutation relations (||). It follows that the two-body wave 
function, defined as in (B), satisfies 



in—ip2{t;x 1 ,x 2 j 



2m 



2 



+ 



i\5(xi - x 2 )(- 







+ 
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dx\ dx\ ' " Z/ V9xi ' 9x2 



^ 2 (t;xi,x 2 ) (58) 



Owing to time- and space-translation invariance, one can separate the time and center-of- 
mass coordinates 



Mt;xi,x 2 ) = e -*(^) e ^)(^)/2 n(xi _ X2) . 



(59) 
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The wave function for relative motion satisfies 



*2 Q2 p2 p 

+ h— X5(x) ) u(x) = Eu(x) . (60) 



V m dx 2 Am 2m 

The presence of the total momentum P in the (^-function potential for relative motion 
vividly demonstrates the absence of Galileo invariance. Provided 

P 



m 



> (61) 



Eq. ( |60| ) possesses a bound state solution with energy 

Since P/m may be identified with (a multiple of) the classical velocity v, we recognize the 
condition (|6l|) as the same as (|l8|) , which guarantees soliton binding. So we expect that 
there is a relation between the classical soliton and quantum bound states. 

Justification for this can be see from the following argument. If we write (|6^ ) as 

p2 

E = — tit~ (63) 

9 A/f 

iIrJ quantum 



we find that M quantum is given by 



2tti 

-^quantum — ■ (^4) 

1 4" 



Next we conjecture that the classical formula for M, Eq. (47), should be modified in 



semiclassical quantization in the same way as in the conventional nonlinear Schrodinger 
equation: viz. as in (|^) iV 3 is diminished by iV.Jl]] Thus we replace (f47| ) by 

M semiclassica i = mN+ (iV 3 - N) . (65) 

For N = 2, this gives 

Msemiclassicai = 2m(l + — ) . (66) 



4 

which agrees with ( |64| ) at weak coupling, i.e., small A 2 . Although explicit solutions to 
the iV-body quantum Schrodinger equation are not known for N > 2, one may establish 
perturbatively in A that (|65| ) is consistent with quantum bound states, at weak coupling. 

The two-body problem inherits from the full iV-body field theory the latter's symme- 
tries and constants of motion. In particular, the total momentum 



P2=Pi+ P2 (67) 



and the total energy 



H 2 = ^—(Pi+pI) ~ 7T~(Pl +P2)^(x 1 - x 2 ) 
2m 2m 

= ^-P 2 + ^-Pl ~ ^-P2\5{x) (68) 
2m 4m 2m 



(p=i(pi-p 2 ), x = x 1 -x 2 ) 



9 



together with the dilaton charge 

D 2 = tH 2 - l(xipi + pixi + x 2 p 2 + P2X2) 

= tH 2 -\ (XP 2 + P 2 X + xp + px) (69) 
(X = \{ xi + x 2 )) 

are time-independent and close on the algebra || 

[H 2 ,P 2 ] = (70) 
[H 2 ,D 2 ] = ihH 2 (71) 
[P2,D 2 ] = ih\P 2 . (72) 

It is interesting that in the present problem scale invariance does not prevent the forma- 
tion of a bound state. That is because the absence of the Galileo invariance allows the 
bound state energy to depend on the total momentum, which is not quantized, and as 
a consequence the bound state energy ( |62| ) lies in the continuum, as is required by scale 
invariance. 

7 Conclusion 

While our chiral, nonlinear Schrodinger equation apparently cannot be completely inte- 
grated, it should be possible, at least numerically, to study the iV-soliton solution. From 
this one could extract the nature of the solitons' mutual scattering and possible dissipa- 
tion. It should then be most interesting to examine the same processes in the iV-body 
quantum Schrodinger equation. Also the procedure for quantizing the classical solutions 
could be further developed. While the analytic intractability may be a daunting obstacle 
for arbitrary N, one should at least unravel the N = 2 case. 
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